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REFLECTED  SHOCK  INITIATION  OF  A  CHEMICAL  PEACTION 

ABSTRACT 

The  Lax  finite  difference  method  is  used  to  compute  the 
hydrodynamic  flow  which  results  from  the  initiation  of  a  chemical 
reaction  by  a  reflected  shock  wave.  The  chemical  reaction  is  assumed 
to  be  irreversible  and  of  first  order;  initial  conditions  are  chosen 
such  that  negligible  reaction  occurs  behind  the  incident  shock  front. 
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LIST  OF  SYMBOLS 


a#  a  c*/V*,  Lagrangian  sound  speed 

a*V* 

a  =  a*/a*  c  -^T 

o 

c*  sound  speed 

c  *  c#/c* 

E*  specific  total  energy  (internal  plus  kinetic) 

E  *  B*/(c£)2 


&  characteristic  length 

m 

mfl.fi s  c oor dinate 

P* 

pressure 

P 

»  p*/p* 

T* 

temperature 

T 

-  T*/TJ 

t* 

time 

t 

»  c*  t*/L* 

0  0 

u* 

flow  velocity 

u 

«  u*/c* 

0 

V* 

specific  volume 

V 

-v*A0 

X* 

distance 

X 

-**/!$ 
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LIST  OF  SYMBOLS  (Cont'd) 
specific  heat  ratio  of  species  A 


y ^  specific  heat  ratio  of  species  B 

7)  mass  fraction  of  A 

Superscript 

*  denotes  dimensional  quantity 

Subscript 

o  denotes  initial  conditions  in  the  unshocked  species  A. 
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INTRODUCTION 


The  study  of  chemical  kinetics  "by  reflected  shock  techniques  has 

1-5* 

been  the  subject  of  many  recent  investigations  .  The  major  advantages 
of  the  reflected  shock  technique  are  that  higher  temperatures  are 
attainable  than  with  incident  shock  techniques  and,  secondly,  the  gas 
behind  a  normally  reflected  vave  is  nearly  stationary,  thus  the 
reaction  is  more  readily  observed. 

In  kinetic  stuii?'*,  'vily  error  in  the  temperature  is  liable  to  cause 
serious  errors  in  the  determination  of  the  reaction  rates;  thus,  it  is 
necessary  to  be  able  to  infer  temperature  precisely  from  measured  shock 
parameters.  Numerous  investigators  have  inferred  shock  temperatures 
behind  reflected  waves  from  measurements  of  pressure,  density,  etc.,  and 
have  found  that  temperatures  may  be  lower  than  theoretical  by  about  2 
percent^  .  Strehlow  and  Case"3  and  Rudinger  ,  however,  have  found  that 
temperatures  behind  the  reflected  wave  may  be  slightly  higher  than  the 
values  computed  from  ideal,  steady-state  shock  tube  theory.  Johnson  and 
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Britton'  have  demonstrated  the  existence  of  lower  reaction  rates  --  which 

imply  lower  temperatures  --  behind  the  reflected  shock  wave,  whereas, 
q 

Fishburne  et  al  have  shown  the  existence  of  slightly  higher  reaction 
rates  behind  the  reflected  wave. 

4 

Strehlow  and  Cohen  have  demonstrated  the  usefulness  of  reflected 
shock  techniques  in  the  study  of  Initiation  of  detonations,  ftiey  have 
observed  the  reaction  wave  behind  the  reflected  shock  and  have  noted 
that  under  certain  conditions  the  reaction  vave  either  developed  into 
a  detonation  wave  before  interacting  with  the  reflected  wave,  or  over¬ 
took  the  reflected  wave  first  and  then  developed  into  a  detonation  wave, 
or,  finally,  the  reaction  merely  accelerated  the  reflected  wave. 


Superscript  nmbere  denote  references  which  may  be  found  on  page  31. 
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EQUATIONS 


Assume  that  ve  have  a  semi -infinite  tube  extending  from  x*  **  -»  to 

x  ■  L*  and  that  the  tube  contains  a  gaseous  species  A  which  is  capable 
o 

of  undergoing  the  irreversible,  exothermic  reaction  A  -•  B.  Furthermore, 

assume  that  the  initial  pressure  distribution  to  be  the  piecewise 

uniform  state,  p£  for  x*  >  0  and  >  p£  for  x*  *  0.  The  flow  velocity 

u*  for  x*  >  0  is  assumed  to  be  zero,  the  temperature  to  be  T  and 
o  o 

specific  volume  to  be  V*.  Initial  conditions  for  x*  £  0  are  given  by 

the  ideal  gas  Hugonlot  conditions.  In  our  sample  calculations  we  shall 

assume  that  p*/p*  *  4.5j  thus  the  incident  shock  speed,  in  units  of  the 

ambient  sound  speed  c*  of  species  A  at  p£,  T£,  and  of  specific  heat 

ratio,  7  ,  of  7/5,  is  2.  The  parameters  in  the  assumed  reaction  rate 

equation  are  chosen  such  that  negligible  reaction  occurs  behind  the 

incident  shock  front.  Under  these  conditions,  the  incident  wave  will 

reach  the  closed  end  of  the  tube,  x*  »  L*,  at  about  L*/(2c*)  time 

o  o  o 

units.  The  shock  reflection  results  in  an  increase  of  temperature 
which  causes  further  reaction  to  occurj  this  increased  reaction  results 
in  the  formation  of  a  compressive  wave  which  eventually  overtakes  the 
reflected  wave.  This  process  will  be  evident  in  the  numerical  solution 
of  the  partial  differential  equations  which  describe  the  reactive  flow. 
See,  for  example,  Figure  12. 


The  equations  describing  the  non-steady,  non -viscous,  one- 
dimensional  flow  of  a  reacting  fluid  in  which  the  irreversible,  first 
order  chemical  reaction  A  -*  B  takes  place  may  be  written  in  the 
Lagrangian  lorn: 


dV  du 

5t  "SS 

du  1  dp 

dt  "  ‘  r0  Sm 


(1) 

(2) 


(5) 
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-  vt]  exp  { -E#/T} 


where  T)  is  the  mass  fraction  of  species  A,  v  is  a  frequency  factor, 

JL  2 

and  Et  is  the  activation  energy,  in  units  of  (c*)  .  The  remaining 
symbols  have  been  defined  previously,  (page  7  ). 

In  addition  to  these  equations,  we  have  the  equations  of  state 


P  •  p(I,V,1|), 


1  2 


where  I  «  E  -  ^  u  is  the  specific  internal  energy,  and 

T  =  T(p,V) .  (7) 

In  the  following,  we  assume  that  both  species  A  and  B  are 
polytropic  gases  with  their  respective  internal  energy  functions  given 
by 

I.  *»  ■  ■  t . +  if  and 

A  W  A 


B  VV3 


+  r* 

T  -D  > 


where  y  ,  ,  if  and  if  are  constants.  Furthermore,  we  assume  that  at 

o'  1'  A  n 

any  given  value  of  p,  V,  and  T  the  mixture  of  A  and  B  obeys  the  ideal 


gas  law 


pV  -  T 


and  that  the  internal  energy  of  the  mixture  of  A  and  B  is  given  by 

I  *  ♦  (1-71)  . 

With  these  assumptions,  equation  (6)  becomes 

f  r  4  ATI  70(70-l)(7i-l)  i 

P  "  i1  “  L^A  *  <^4]}  (7,-1  W  ♦  (l-^(70-l  (VJ 

and  T  is  given  by  equation  8. 
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difference  equations 


The  partial  differential  equations  (l)  -  (5)  must  be  replaced  by 
a  system  of  difference  equations  in  order  to  obtain  a  numerical 
solution.  The  equations  (l)  -  (5)  are  all  of  the  form 


tf  +  ijs  +  B 

at  am 


0 


where  F  and  B  are  functions  of  m,  t  and  y.  Such  equations  are  called 
conservation  laws  and  are  of  the  form  used  in  the  difference  method 
proposed  by  Lax10  in  1954 .  In  the  Lax  method  time  derivatives  are 
replaced  by 


at 


yd  -  i  (yd +  yd"d-\ 

yn  2  vyn-l  yn+l' 


At 


and  space  derivatives  F^  by 


BF 

Si* 


2fm 


where  superscript  J  denotes  the  value  of  the  variable  at  time, 
t  =  td”1  +  Atd  and  subscript  n  refers  to  the  spatial  position  m  **  nAm. 
The  spatial  grid  size  Am  is  fixed  whereas  the  time-step  depends  upon  J. 
Letting  F  **  ^  (F^"^  +  F^_^)  the  difference  equations  approximating 
equations  (l)-(5)  become 


rd 

n 

J 

n 


Ed  *  E 
n 


V  +  ~  (u d"|  -  uJ'j) 

2  '  n+1  n-1 

(10) 

3  -  ^ 

(n) 

'  [(pu)n+l  '  (pu)nll] 

(12) 

|  -  vAtd  If  exp  (E^/pv) 

(U) 
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*n  ‘  *  +  ¥  R  +  S]  (14) 

i  "  En  '  I  ("n>2  (15) 


PJ  «  1- 
11  VJ 

n 


K  -  GJ4 +  (16) 


M  i  ♦  (v1)]1 

where  X^  =  At^/Anu 

The  ratio  X^  =  At^/Am  is  calculated  from  the  stability  condition  that 


max 

n 


i-iY 
*  aJ  • 
n 


for  0  <  b  *  1 


and  where  a,  the  Lagrangian  sound  speed,  is  given  by 

1 1/2 


a 


J 


71(70-!)  +  (V’vX  Pn 


J 


"oK-^n +  v4  vi 


INITIAL  AND  BOUNDARY  CONDITIONS 

In  terms  of  the  non-dimensional  variables  previously  defined,  the 
mass  of  gas  A  contained  in  the  region  0  <  X  s  1  is  equal  to  one.  This 
region  is  divided  into  N  equal  mass  zones,  i.e., 


hence,  for  this  section  the 
the  same  numerical  value  as 
-»  <  X  *  0  is  replaced  by  a 
The  initial  conditions  are: 


Lagrangian  mass  coordinate  has,  initially, 
the  Eulerian  coordinate.  The  section 
finite  section  containing  M  mass  zones. 
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RESULTS 


The  results  of  the  numerical  computation  of  equations  (10)  -  (15) 
are  summarized  in  Figures  1-7.  Figure  1  is  a  space-time  plot  of 
selected  particle  paths  and  shows  the  general  features  of  the  flow. 

In  particular,  the  incident  and  reflected  shocks  and  the  reaction  wave 
are  clearly  evident.  After  interaction  of  the  reaction  wave  and 
reflected  wave  occurs,  a  double  wave  appears.  (The  double  wave  is 
more  evident  in  Figure  12).  Each  of  the  two  waves  appears  to  be  travel 
ing  at  nearly  constant  velocity.  The  secondary  wave,  receding  from 
the  first,  is  a  shock  wave  traveling  in  gas  B.  The  first  wave  appears 
to  have  adl  the  characteristics  of  a  detonation  as  the  values  of  the 
pressure,  etc,  correspond  to  those  expected  from  a  detonation 
transition  of  species  A  into  species  B,  under  the  chosen  initial 
conditions.  The  observed  double  wave  structure  is  similar  to  that 
observed  by  Cher'  and  Kistiakowsky11  in  their  photographic  studies  of 
detonation  of  certain  hydrocarbon-oxygen  mixtures .  Cher  and 
Kistiakowsky  concluded  that  "the  secondary  wave  is  due  to  entropy 
increase  in  the  rarefaction  wave  caused  by  a  spontaneous  reaction”. 
Analysis  of  the  wave  patterns  observed  in  these  computations  will  be 
reported  in  subsequent  reports,  and  it  is  expected  that  further 
details  concerning  the  double  wave  structure  will  be  reported  at  that 
time. 

A  typical  pressure -distance  curve,  corresponding  to  a  fixed  time, 
is  shown  in  Figure  2.  This  pressure  profile  is  typical  of  the  pressure 
wave  after  the  occurrence  of  interaction  between  the  reflected  and 
reacting  wave.  Negligible  reaction  has  occurred  ahead  of  the  primary 
wave  (corresponds  to  peak  pressure)  and  the  conditions  in  this  region 
are,  approximately,  the  same  as  behind  the  initial  incident  shock 
front,  i.e.,  p  -  4.5,  u  ■  1.25,  etc.  An  ideal  Chapman -Jouguet 
detonation  wave  traveling  into  species  A  at  these  conditions,  and 
instantaneously  converting  A  into  B  gives 
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t 


PCJ  «  84.70  Vgj-  0.2237 

-  -1.69^  D  -  -6,047 

vbere  D  denotes  the  speed  of  the.  wave.  The  computed  values, 
corresponding  to  the  peak  pressure  of  Figure  2,  are 

p  -  84.62  V  -  0.2239 

u  *  -I.692  D  •  -6.045 

The  tail  of  the  pressure  profile  of  Figure  2,  i.e. ,  the  nearly 

constant  pressure  region  adjacent  to  the  vail,  corresponds  to  the 

region  vhere  the  particle  velocity  Is  approximately  aero.  If  ve  assume 

that  this  region  is  connected  to  the  detonation  state  by  means  of  a 

rarefaction  vave,  the  pressure  ratio  at  the  vail  vould  be  48.02  as  • 

compared  with  the  computed  pressure  ratio  of  48.17.  The  vail  pressure 

at  J  »  1049,  300  time -steps  later,  Is  48.02  and  the  quiescent  region 

occupies  a  corresponding  greater  region  of  the  tube.  Comparison  of 

12 

these  results  vith  0.  I.  Taylor*  s  similarity  solution  vlll  be 
reported  In  subsequent  reports. 

If  ve  assume  that  the  region  bounded  by  the  vail,  the  reflected 
shock  and  the  reaction  vave  is  constant  and,  In  particular,  assume 
the  values  of  pressure  and  velocity  are  those  given  by  the  Ideal  gas 
shock  equations1^,  then  for  an  Incident  shock  pressure  of  4.5,  the 
pressure  and  particle  velocity  In  the  reflected  zone  are  p  ■  15  and 
u  ■  0.  The  Chapman -Jouguet  detonation  vave  values  vould  be 

PCJ.  198.1  Vqj  ■  0.1004 
Upj  -  -2.944  D  -  -7.403  . 

For  comparison.  Figure  3  shows  the  pressure,  at  constant  time,  Just 
prior  to  the  Interaction  of  reflected  and  reaction  vaves.  The 
pressure,  corresponding  to  the  peak,  is  199*7  and  the  corresponding 
values  of  u  and  V  are 


/ 
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U  a  -3.059> 

V  a  0.1003, 

and  the  reaction  is  traveling  at  speed  D  »  -7.504. 

Conditions  behind  the  reflected  wave,  prior  to  the  observed 
development  of  the  reaction,  are  not  constant  as  can  be  seen  in 
Figures  4  and  5  which  are,  respectively,  plots  of  the  wall  pressure 
and  wall  temperature  versus  time. 

The  reaction  is  complete,  i.e.,  T)  =  0,  at  the  time  corresponding 
to  the  maximum  pressure  of  Figure  3.  The  mass  fraction,  T),  and  the 
specific  volume  ratio,  at  X  =  1,  are  shown,  respectively,  in  Figures  6  and  7. 

If  v  in  equation  12  is  set  equal  to  zero  and  the  remaining 
parameters  and  conditions  are  unchanged,  the  problem  considered  is 
merely  that  of  normal  shock  reflection  in  an  ideal  gas.  For  purposes 
of  comparison  and  preliminary  tests  of  accuracy  v  was  set  equal  to 
zero  and  the  corresponding  flow  computed.  In  Figure  8  the  particle 
paths  are  plotted.  The  results  are  those  expected,  i.e.,  the  particle 
velocity  is  (nearly)  zero  behind  the  reflected  wave  and  the  pressure 
Jump  is  that  given  by  the  well-known  shock  reflection  formulas1^.  The 
pressure  as  a  function  of  Eulerian  distance  is  plotted  for  various 
grid  sizes  in  Figure  9>  and  it  is  seen  that  the  shock  front  becomes 
steeper  as  grid  size  decreases. 

In  Figures  10  and  11,  percent  relative  errors  In  V  and  X  are 
shown  as  functions  of  grid  size.  The  errors  in  p  and  u  are  negligible, 
at  these  chosen  points,  for  all  grid  sizes. 

The  effect  of  increasing  the  activation  energy  on  the  resulting 
reacting  flow  is  illustrated  in  Fig.  12.  The  activation  energy  E*  was 
increased  from  25  to  28.6.  The  values  of  the  remaining  parameters  are 
the  same  as  those  used  to  obtain  the  results  illustrated  in  Fig.  1. 
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Fig.  2  Pressure  Ratio  vs  Scaled  Distance 
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